Parametric resonance for antineutrino conversions using LSND best-fit results with a 

3+1 flavor scheme 
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An analytical solution to a parametric resonance eflFect for antineutrinos in a 3+1 flavor (ac- 
tive+sterile) scheme using multiple non-adiabatic density shifts is presented. We derive the con- 
ditions for a full flavor conversion for antineutrino oscillations Va — > Fs (q = e, fi, r) under the 
assumption that LSND best-fits for the mixing parameters are valid in a short-baseline accelerator 
experiment. We show that the parametric resonance effect can be exploited to increase the effective 
antineutrino oscillation length by a factor of 10-40, thus sustaining a high oscillation probability for 
a much longer period of time than in the vacuum scenario. We propose a realistic experimental setup 
that could probe for this effect which leaves a signature in terms of a specific oscillation probability 
profile. Moreover, since the parametric resonance effect is valid in any 2 or 1-1-1 flavor approximation, 
our results could be suggestive for future short-baseline accelerator neutrino detection experiments. 

PACS numbers; 13.15.-|-g, 14.60.Lm, 14.60.Pq 



I. INTRODUCTION 



Ever since Pontecorvo gave birth to the concept of 
neutrino oscillations in 1958 jL[, the model has evolved 
significantly over the years [2 Q , allowing us to have a 
comprehensive understanding of neutrino behaviour to- 
day. Oscillations take place in both vacuum and matter. 
The presence of a medium gives rise to matter potentials 
felt by the neutrinos propagating through it 0, IE IE 
which in short alters the oscillation probability compared 
to propagation in vacuum. This effect can have dra- 
matic consequences in resonance situations. For a non- 
adiabatic density shift, the exact analytical solution for 
the oscillation probability is available in e.g. Ref. 
while various uses of a single and multiple shifts have 
been studied in e.g. Refs. |EE,0|- The idea of using 
multiple unit-step density shifts to obtain a full flavor 
conversion for neutrinos was first studied in Ref. 12J. 

Although the scenario of neutrinos propagating 
through several layers of matter with constant densities 
certainly is not new, this paper contributes to this par- 
ticular topic through the study of resonance conditions 
for antineutrino flavor conversions in a 3-1-1 flavor sce- 
nario with LSND best-fit parameters. One could very 
well ask for a justification concerning the choice of a 3+1 
flavor scenario as there certainly has been voices in favor 
of for instance a 3+2 solution [lE|. Adding one sterile 
neutrino flavor represents the simplest extension of 
the Standard model, containing an active neutrino fla- 
vor for each lepton generation e, /i, r. The need for such 
an extension comes from the discrepancy between solar 
and atmos phe ric neutrino detection experimental data 

EE EE El 

llSl and the short-baseline accelerator LSND 
|l9| results if one insists on using a 3+0 scheme. Consis- 
tency can be regained by adding a sterile neutrino flavor, 



thus obtaining three independent square mass differences 
Am^. The opportunity to signiflcantely influence the 
oscillation probability and oscillation length arises in a 
3+1 flavor scenario under valid approximations concern- 
ing the square mass differences, as we will show. This 
the main motivation for selecting the 3+1 flavor scheme. 

This paper is organized as follows. In Sec. ^Iwe es- 
tablish the 1+1 flavor approximation for LSND best-fit 
results, while in Sec. Illll the parametric resonance effect 
is treated for this scheme. The opportunities that arise 
from the parametric resonance effect with regard to ex- 
ercising control over the neutrino oscillation probability 
and oscillation length are discussed in Sec. IIVI where we 
also show that this effect would be observable in a real- 
istic experimental setup. Concluding remarks are given 
in Sec. |V| 



II. OSCILLATIONS IN THE 3+1 FLAVOR 
SCHEME 



The neutrino flavor states {va) (ct — e, /x, t, s) are linear 
superpositions of the mass eigenstates jt'i) {i = 1,2,3,4) 
and vice versa. The mixing between these two sets of 
states explicitely reads 



(1) 



1=1 



where J/ is a unitary 4x4 neutrino mixing matrix which 
we parametrize in a similar manner as Ref. |T^ . 



U = U3iU24U23UliUuUi2. 



(2) 
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The matrices Uij represent mixing between pairs of mass 
eigenstates and have the generic form of a rotation with 
an angle 0ij in 2-dimensional i — j space. For example. 
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one has 
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(3) 



our analysis, this translates into /S.mf^^^^ dominating the 
mass scale, thus also dominating the oscillation arena. 
This justifies an approximative model of two flavors for 
the LSND short-baseline accelerator experiment. We are 
immediately able to re-write Eq. (jsj as the familiar for- 
mula 



where = cos Oij, Sij = sin 6ij. In writing down [/, we 
have intentionally left out the CP-violating phase ^cpj 
assuming it is sufhciently small to be neglected, and Ma- 
jorana mass phases, which are irrelevant with respect to 
oscillations. The time dependence of the mass eigen- 
states is governed by |i/i(x, i)) = e~'-^'*|t'i(x, 0)), where 

ii^/2E, with the definition 



E 



IpI = E. Consequently, we operate with neutrino mo- 
menta p much greater than their masses m^, which is 
an excellent approximation, and assume equal momen- 
tum among the mass eigenstates. Furthermore, we adapt 
the mass hierarchy scheme in Fig. ^with the definitions 
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(6) 



Having gone from 3-|-l to 1+1 flavors, an interest- 
in^ opportunity arises. The parametric resonance effect 
p2] is now directly applicable in the approximation that 
AmgoL f-iid ^"^ATM negligible compared to AmLgND- 
Due to the large value of Aml^-^-^ , the corresponding os- 
cillation length is relatively short; of order 10-100 m for 
neutrinos with E ~ 10-100 MeV. It is important to realize 
that since Aml^,^^ = Am^^ ~ Aml^ — Aml^, the fol- 
lowing analysis is applicable to the Va — > (a = e, /i, r) 
channel. The consequences of an enhanced rate of flavor 
conversion from Va Vs could be tested in a disappear- 
ance experiment, where a lack of Va signatures would 
indicate the occurence of such oscillations. 



m2 
mi 



FIG. 1: Assumed mass hierachy scheme for neutrino mass 
eigenstates \vi{-x.,t)). 



III. PARAMETRIC RESONANCE CONDITION 

We now explicitely derive the parametric resonance 
conditions for a full antineutrino flavor conversion in the 
Va Vs channel. Consider an experimental setup as in 
Fig. 12 Since oscillations are altered in the presence of 



The neutrino oscillation probability is calculated 
in standard quantum mechanical fashion Pa-^p — 
\{i^p{x,t)\i^a{x,t))\'^, producing 
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Neutrino production 
point t = 



XI 



d AL 



AL 



AL 



with Jafjij = UpiU*^Up^Uja. The orthonormality condi- 
tion (i/i(x, i)|i'j(x, t)) = 6ij is assumed to hold. 

Considering the mass hierarchy scheme of Fig. ^ it is 
seen that Eq. Q can be significantely simplified. The 
best-fit results from the LSND experiment states that 
AttIlsno ~ 1.2 eV^ 0, such that it is a very good ap- 
proximation to set AttIlsnd » {Aml^^, Am\^^} since 
the best-fit equivalents for these are 7.1 x 10~^ eV^ and 
2.5 X 10-3 eV^, respectively ol El El. Upon this 
approximation, Eq. turns into (see e.g. Ref. E3l) 



Pa^f3 = Safj - 4:Ua4Uf34^af3 - UaAU^i] siu' 
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(5) 



Eq. Q is nothing but the ordinary oscillation formula 
in the 2-flavor scenario in disguise. One can understand 
this effective reduction from 3-1-1 to 1-1-1 flavors qualita- 
tively by what is referred to as the "quasi two neutrino 
approximation" or " one mass scale dominance" |2^ . For 



FIG. 2; Using multiple non-adiabatic density shifts to obtain 
a flavor conversion Va Vs. Material slabs of density p 
and length d are placed at regular intervals in the neutrino 
trajectory. 

matter, we must take neutrino forward scattering reac- 
tions with particles into account. A list of such matter 
potentials can be found in Ref. [23 ■ For oscillations 
between active flavors, it is a known result that only 
the charged-current (CC) contribution from scattering 
on electrons is of significance with respect to oscillations 
since the neutral-current (NC) contribution just gives an 
overall phase which is irrelevant. In the present case, 
we consider oscillations between one active and one ster- 
ile flavor. The equation of motion for the antineutrinos 
thus reads 



dt 



= (Ho + Hi) 



(7) 
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where Hj contains the forward scattering potentials and 
is given by 







(8) 



The CC and NC potentials read = V2Gj.A^„ and 
= —GpNn/ y/2 0, where Gp is the Fermi constant, 
while Na and Nn are the number densities of leptons 
with flavor a and neutrons in the medium, respectively. 
The minus sign in Eq. (jS]) comes from the fact that the 
potential sign for antineutrinos is reversed compared to 
neutrinos |2l| . Accordingly, antineutrinos under the in- 
fluence of a potential must be described by a different set 
of mixing parameters in the medium than in vacuum. In 
short, the transformation from the vacuum scenario to 
the presence of matter is obtained by substituting the 
vacuum parameters with their matter equivalents, i.e. 
6 dm and Am^ — + AAf^, for the case of a constant, 
isotropic electron density. In the present case, these read 
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where the notation 



(9) 
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If a non-adiabatic density shift occurs at t = to, the 
neutrino propagator that takes the neutrino state vec- 
tor from t — Q to t — to reads 



IC{t) = C/diag(e 



-iEat 



-lEit 
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(12) 
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Here, [/ is a 2 x 2 mixing matrix corresponding to a 
rotation with an angle Oas in the mass eigenstate-plane 
for Va and 77^, such that 



(13) 



As pointed out in Ref. , the flavor eigenstates must be 
continuous in the transition between vacuum and a mas- 
sive medium even if the density change is extremely non- 
adiabatic. Demanding that the flavor eigenstates are con- 
tinuous across the density shift, we write |4'(x, to))™^^^ = 
C/,^^|^(x, to)), where the subscript m denotes our transi- 
tion to a massive medium. The propagation through the 
medium is described by replacing the vacuum quantities 
with their medium equivalents, such that the flavor state 
vector in the medium is 

|*(x,t,io)> =/C™(t-io)'»C(to)|^'(x,0)), t>to. (14) 

Here, JCm{t — to) is defined in a similar fashion as the 
vacuum propagator, namely 



has been introduced. For short-baseline accelerator ex- 
periments, one can effectively set Na = for a S {^, t} 
such that only the NC potential contributes. For a — e, 
however, one must include Ne in Eq. H1U|) . 

Consider now the flavor state vector |'I'(x, i)) describ- 
ing the antincutrino state. Our initial condition is that 
|^'(x, 0)) = [1 0]'^, using a basis where 

ii7„(x,o)) = [ior, ii7,(x,o)) = [oir. (11) 



ICrnit) = ?7mdiag(i 



u- 



(15) 



Eq. thus gives us the opportunity of using the 

IC{t) and ICm{t) propagators to "rotate" the original fla- 
vor eigenstate Fq into a Vs by an appropriate choice of 
parameters. For a setup such as in Fig. [3 we obtain 



|«'(x, Lo, AL, n, d)) = ICroALo, AL, n, d)\^{x, 0)), 



(16) 



(n-l) 



/CTOT(io,Ai,7i,d) -/C„,(d) X ^/C(AL)/C„(d)J /C(Lo), n= 1,2,3,.. (17) 

Now, Eq. (|17|l can be written 

|*(x,Lo,Ai,ri,d)) = [(7™D™(d)C/-iC/i?(AL)t/-i]("-^)t/™i^™(d)(7-i;72?(Lo)C/"'|*(x,0)), (18) 

where D{x) = diag(l, e~"^™LSND^/2£^) and D,n{x) = diag(l, e""^*^^^/^^). To obtain the diagonal matrices D and D^, 
we have extracted common phase factors of the type e"^ which are irrelevant when pursuing the oscillation probability. 
The resonance |5'(x, Lo, AL, n, d)) = [0 1]'^ is obtained when we set Lq = AL and choose the phases such that 



Amls^^AL/2E = n, AM^d/2E ^ tt. 



As a consequence, D{AL) = D„i{d) = diag(l,— 1), and one obtains 
[UmD„,{D)U,-^][UDiAL)U-^] = 



cos 2{9n 
-sin 2 (61 



U) sin2( 
6*04) cos2( 
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(19) 



(20) 
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FIG. 3; Oscillation probability Pose for the setup in Fig. |21as 
a function of the number of iterations n. The parametrically 
enhanced oscillation probability reaches unity at n = 7r/4AS 
for i\E = 0. 



where is the Pauh matrix 



-i 

1 



(21) 



Eq. (^SJ) then becomes |*(x, Lq, AL, n, d)) 

g2m(e„-e„4)^''|^(x,0)), written out as 



|4'(x, Lo,AL,n,d)) 



cos(2nA6') 
- sin(2nA6') 



(22) 



with the definition A9 = {6m — dai)- The resonance 
condition for the number of iterations n is then 



2nA6l = ±7r/2. 



(23) 



If this condition is met, the probabihty for a flavor con- 
version reads Pas = sin^(2nA0) = 1. It is important to 
reahze that this resonance cruciaUy depends on choosing 
the phases properly. To see this, recaU that d and AL 
must be chosen to satisfy Eq. H19|l . If d and L instead 
are shghtly perturbed to satisfy 



4,„AL/2(i? + A£;)=7r 
AM^d/2{E + AE) =7r. 



the conversion probability never reaches unity and is 
displaced. This is shown in Fig. |31 In general, the width 
of the resonance depends on A9. For AE/E ^ Ad, it 
is seen from the graph that one still gets reasonably 
close to the resonance. Larger values of AE destroys 
the resonance, and gives a flat curve for the oscillation 
probability. 



IV. DISCUSSION OF RESULTS 

We showed above that it is, in theory, possible to ob- 
tain a full conversion of Va — > (a = e, /i, t) in a short- 
baseline accelerator experiment with LSND best-fits for 
the neutrino mixing variables given a resonance set of ex- 
perimental parameters. While the required material slab 
lengths are manageable, there is also the issue of mate- 
rial density. When restricted to realistic densities in the 
range p — [1, 11] g/cm'^, one obtains A6 ~ 10~^ for the 
LSND best-fit ^LSND — 0.019, resulting in an extremely 
narrow resonance. The number of iterations required to 
reach the resonance of Pas will thus be very high unless 
there is a way of obtaining a large A9 in the medium. 

Although the idea of obtaining the full resonance for 
antineutrino conversions is unrealistic for material den- 
sities available in an experimental setup, an interesting 
opportunity arises from our analysis. For vacuum and 
matter oscillations in a single material, the oscillation 
length for an ii^ = 50 MeV neutrino using LSND pa- 
rameters is approximately 50 m. We now show that one 
is able to prolong the effective oscillation length to sev- 
eral hundred meters by using the parametric resonance 
effect in an experimental setup with material slabs of re- 
alistic density and length. In this context, the effective 
oscillation length refers to the distance upon which the 
resulting oscillation probability pattern completes a full 
period. Note that this is neither the vacuum oscillation 
length ^vac nor its matter equivalent Im, which of course 
remain unaltered. Accordingly, the parametric resonance 
effect would constitute a way of manipulating the neu- 
trinos to maintain a large oscillation probability for a 
longer period of time. Such an effect could be probed 
for with aim to investigate the validity of the 1-1-1 flavor 
approximation for LSND parameters, and could also be 
suggestive in future neutrino detection experiments using 
short-baseline accelerators. 

Consider Fig. ^for numerical simulations of the para- 
metrically enhanced oscillation probability for Ve <-> Vs 
oscillations. Several plots are included for different 
choices of d and L, i.e. deviations AE in Eq. H24|l . 
Note that this effect does not occur only within a narrow 
window of deviations AE/E ^ AO which was the condi- 
tion for the full resonance. Here, deviations as large as 
20% would still cause the oscillation length to be several 
times larger than in the vacuum scenario. For densities 
p G [1,11] g/cm^, the numerical results in Fig. ^ are 
essentially unchanged. As a consequence, the results are 



(24) also equally valid for 



and 



oscillations 



since the effect of including in Eq. I|10(l is negligible, 
i.e. A9 remains practically unaltered. 

Although the maximum oscillation probability in Fig. 
^ is the same as for vacuum oscillations, i.e. Pmax = 
sin^ 20LSND — 0.0015, it is clear that the effective oscilla- 
tion length increases by a factor of 10-40 when using the 
setup of Fig. 12 with realistic material slab densities. In 
comparison, neutrino propagation through a single ma- 
terial slab of density p = 3 g/cm'^ has virtually no ef- 
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feet on the oscillation length due to the smallness of A6'. 
In order to probe for this parametric resonance effect, 
one would consequently need an experimental setup with 
consecutive material slabs of density p ~ 3 g/cm'^ and an 
effective neutrino trajectory of less than 1 km. We stress 
the fact that the idea of substantially increasing the oscil- 
lation length by means of the parametric resonance effect 
is experimentally realistic and constitutes an interesting 
meta-effect that could yield information concerning the 
validity of a 2 or 1+1 flavor approximation in a neutrino 
factory. Thus, the results presented in this paper could 
prove useful in light of future short-baseline accelerator 
neutrino detection experiments. 

Finally, we comment upon potential obstacles and 
competing effects that are vital to take into account 
when considering this type of parametric resonance ef- 
fect. For instance, in a disappearance experiment where 
a deficit of Da signatures is measured, it is important to 
estimate the magnitude of Va conversions where 

/3 ^ s in order to be certain of that one is truly ob- 
serving the present parametric resonance effect. Con- 
sider for instance production of at a neutrino source 
in a disappearance experiment where 9^, oscilla- 

tions would be a competing effect with the parametric 
resonance presently discussed. Making use of Eq. Q to 
estimate the conversion of active-active type as opposed 
to active-sterile, it is found that for the established large 
mixing angle (LMA) parameters aforementioned for dis- 
tances L ~ 1 km, Pp^_^p_, = 5.5 X 10^^ when the neutrino 
energy is = 50 MeV. This amounts to roughly one 
third of the maximun oscillation probability in Fig. 0] 
which obviously is a non-negligible contribution to the 
total conversion rate of P^. Discrimination of active- 
active conversions from active-sterile conversions could 
be obtained from e.g. the process -I- p — > e+ -I- n in 
a scintillator detector. We stress that the active-active 
flavor oscillations are not subject to the parametric res- 
onance conditions given our present values of Lo,d, and 
AL (see previous section) and simply mimic the ordinary 
vacuum oscillation probability in spite of the presence of 
slabs. The active-sterile flavor oscillations, however, are 
subject to the interesting manipulation illustrated in Fig. 
0] Moreover, it is clear that the observation of this effect 
demands that the speciflc oscillation probability proflle of 
Fig. ^is identified with enough accuracy to unquestion- 
ably distinguish it from pure vacuum oscillations. This 
would preferably require several detectors placed along 
the trajectory of the antineutrinos instead of one; e.g. 
an intermediate detector at for instance L ~ 0.5 km. 
Clearly, the solid line shown in Fig. ^ offers the most 
promising scenario for the study of the oscillation proba- 
bility proflle since fewer oscillations will have taken place 
than in the case of the dashed and dotted lines at any 
distance. Another issue that also deserves comment con- 
cerns the collimation of the neutrino beam. A large angu- 
lar spreading of the neutrino trajectories could endanger 
the validity of the parametric resonance conditions for a 
setup as shown in Fig. |21 Short-baseline experiments of 




Uot [m] 



FIG. 4: Parametrically enhanced probability for Ve <-» Vs os- 
cillations using i? = 50 MeV LSND antineutrinos propagating 
tlirough material slabs with density p = 3 g/cm^. The oscil- 
lation length is significantely increased compared to vacuum 
oscillations and matter oscillations in a single material slab. 
For T^ct <-> 77s, Q G {p, t} oscillations, these numerical results 
are also valid since remains practically unchanged. With 
a total neutrino trajectory of less than 1 km, a realistic ex- 
perimental setup would be able to probe for this effect. 



type LSND with a distance source-detector of less than, 
say, 100 m, will beneflt from a better collimation than 
distances approaching 1 km. The collimation in the lat- 
ter case would also be worse than neutrinos with ener- 
gies of order ©(GeV) in long-baseline experiments. As 
an estimate for the loss of neutrino flux poor collimation 
can give rise to, thus lowering the number of neutrino 
trajectories originating from the source that arc subject 
to the parametric resonance conditions, consider the fol- 
lowing scenario. Assume that the slabs in Fig. [21 are 
of height 2 m and that the entire neutrino flux is inci- 
dent on the flrst slab located at Lq but spreads out in 
a cone from the source. At a distance of L = 1 km, 
the neutrino flux through the last slab {e.g. a detector) 
would then have been reduced to a fraction / = 
where i?o is the radius of the cross-section of the cone 
at that distance. With the given conditions, one easily 
flnds that _Ro = itan0 where (p = atan(l m/Lo) is the 
angle at which the neutrino beam is spread compared to 
horisontal ground. Thus, / ~ 1/400 and it is seen that a 
considerable fraction of the neutrino flux is removed from 
the trajectories passing through the slabs. Consequently, 
this would constitute a possible limiting factor of the ap- 
plicability of the parametric resonance effect depending 
on the size of the neutrino flux of the source, or at least 
place high demands upon the collimation of the neutrino 
beam. 
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V. SUMMARY 

In summary, we have derived an exact analytical ex- 
pression for parametric resonance for antineutrino flavor 
conversions in a 3+1 flavor scenario for LSND best-fit 
results. In order to exploit the idea of multiple non- 
adiabatic density shifts to obtain an oscillation probabil- 
ity that exceeds the maximum oscillation probability in 
vacuum (matter), it follows from Eq. H22|) that 

n > d^^)/A0 (25) 

must be satisfied. Consequently, the observability and 
application of a parametrically enhanced oscillation prob- 
ability depends on obtaining a large A0, which is im- 
possible to achieve for realistic densities. However, the 
parametric resonance effect offers an interesting possibil- 
ity to increase the effective neutrino oscillation length, 
thus leading to the maximum vacuum oscillation prob- 
ability being sustained over a much larger interval than 
in the scenario of pure vacuum oscillations or even oscil- 
lations in a single medium. We showed that this effect 
would be observable in a realistic experimental setup us- 
ing material slabs of density p ~ 3 g/cm^ and a total 
neutrino trajectory of less than 1 km, giving rise to a 
specific oscillation probability profile where the effective 
oscillation length is manifested as shown in Fig. 01 We 
stress that the maximum oscillation probability is not en- 



hanced over this distance, but rather the effective oscil- 
lation length. With regard to detection, a disappearance 
experiment could search for Vs oscillations in both 

the CC and NC channel, while — > Vs and F,- — > V s 
oscillations would be harder to detect since interactions 
with Fi, i ^ [i^T only occur in the NC channel. Further- 
more, there are experimental complications associated 
with probing of the parametric resonance effect in terms 
of competing active-active flavor oscillations and the col- 
limation of the neutrino beam over a longer distance 
source-detector (1 km) than most short-baseline accel- 
erator experiments. These arguments must be taken into 
account upon considering realistic approaches in terms 
of exploiting the parametric resonance effect for active- 
sterile antineutrino conversions. It is clear that this para- 
metric resonance of the oscillation probability would be 
applicable to any 2 or 1+1 flavor scenario. In light of 
the expected MiniBooNE results, the validity of such an 
approximation could yield an interesting opportunity to 
take advantage of this effect. 
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